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ABSTRACT 

A systematic investigation of dust emission associated with the ionized gas has so far been performed only 
in our Galaxy and for wavelengths longer than 60 fim. Newly available Spitzer data now offer the opportunity 
to carry out a similar analysis in the Large Magellanic Cloud (LMC). 

By cross-correlating Spitzer SAGE (Surveying the Agents of a Galaxy's Evolution) data with the 
ATCA/Parkes HI 21 -cm data, the NANTEN 12 CO (J=l-0) data, and both the SHASSA Ha and the Parkes 
6-cm data, we investigate the physical properties of dust associated with the different phases of the gas (atomic, 
molecular and ionized). In particular, we study the presence and nature of dust from 3.6 to 160 fim and for 
various regimes of the ionized gas, spanning emission measures (EM) from ~ 1 pc cirT 6 (diffuse component) 
to ~ 10 3 pc cirT 6 (HII regions). 

Using a dust emission model, and testing our results with several radiation field spectra, we show that dust 
in the ionized gas is warmer than dust associated with other phases (atomic and molecular). We also find a 
decrease of the polycyclic aromatic hydrocarbons (PAH) relative abundance with respect to big grains (BGs), 
as well as an increase of the near infrared (NIR) continuum. These three results (e.g. warmer temperature, 
decrease of PAH abundance and increase of the NIR continuum) are found consistently for all regimes of the 
ionized gas. On the contrary, the molecular phase appears to provide favorable conditions for the survival of 
PAHs. Furthermore, the very small grain (VSG) relative abundance tends to increase in the ionized phase, 
especially in bright HII regions. Last but not least, our analysis shows that the emissivity of dust associated 
with the ionized gas is lower in the LMC than in our Galaxy, and that this difference is not accounted for by 
the lower metallicity of the LMC. 



1. INTRODUCTION 

Detecting dust emission associated with the diffuse ionized 
phase is a complex task due to two main competing factors. 
First of all, large grains are expected to be selectively de- 
stroyed in the warm ionized medium (WIM) due to the pas- 
sage of multiple shock fronts which would reduce their emis- 
sion. In addition, along most line of sights (LOS), the diffuse 
ionized gas is likely mixed with the other phases of the gas 
(atomic and molecular), and its contribution is not expected 
to be dominant at any latitude. 

Because of these difficulties, contradictory results concern- 
ing the pre sence of dust in the W IM have been reported in the 
literature. Lag ache et al.l (119991) . looking at the high Galac- 
tic latitude residuals of the IR/HI correlation between 100 
and 1000 fim, claim a statistical detection. In particular, 
by adopting an emissivity spectral index of 2, they find an 
emissivity law in the atomic and ionized gas phase of, re- 
spectively, r/N H = 8.7 ± 0.9 10~ 26 (,1/250 fim)- 2 cm" 2 and 
t/N h = 1.0 ± 0.2 10~ 25 (A/250 fim)- 2 cirT 2 , with Tpi = 17.5 
K an d T//+ - 20 K. Fol lowing this preliminary investiga- 
tion, lLagache et al.l (120001) extend their analysis, and carry 
out a decom position of the far in frared (FIR) emission us- 
ing WHAM flRevnolds et alJI 19981) Ho- data as a tracer of the 
ionized medium. The result of this study confirms the early 
detection, with an emissivity law for the ionized phase of 
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t/N h = 1.1 10~ 25 (,1/250 fim)- 2 cm" 2 , for a fixed tempera- 
ture of 17.2 K for b oth the atomi c and i onized phase. At low 
Galactic latitudes, Paladi ni et alj ((2007) derive dust emissiv- 
ities associated with the atomic, molecular and ionized gas 
phases, using an inversion method applied to IRAS (Neuge- 
bauer et al. 1984), DIRBE (Hauser 1993), Archeops (Benoit 
et al. 2002) and WMAP (Bennett et al. 2003) data. Ac- 
cording to their analysis, the average dust temperatures in the 
atomic and molecular phases are comparable, and of order 
of 19.8 and 19.2 K respectively, while dust in the ionized 
phase is significantl y warm er, T#+ 26.7 K. At the same 
time, lOdegard et ail (120071) looking at the DIRBE/WHAM 
correlation in five selected high -latitude Galactic regions find, 
contrary to Lag ache et al.l (120001) . that the Ha-correlated dust 
component is negligible in all cases. 

Despite the observational caveats, assessing the presence of 
dust in the ionized gas phase is very important, given that the 
optical properties of dust grains can vary with the properties 
of the local environment. To date, no cl ear-cut indication of 
such a behaviour has been found though. iParadis et al.l (1201 lb 
perform a statistical study of all the molecular clouds in the 
LMC and do not find any systematic change in dust abun- 
dances nor in dust temperature be tween the mo l ecular and 
the atomic phase. On the contrary, IStepnik et al.l (120031) and 
Paradi s"et al.l (|2009b), by looking at nearby Galactic clouds, 
do report evidence of this kind of variations. In the ionized 
medium, due to the proximity of the ionized gas with UV 
sources such as O and B stars and the cumulative effect of 
shock fronts, changes in the optical dust properties are indeed 
expected. For instance, depletio n of PAHs has been observed 
toward individual HII regions ([Eve rett & Churchweli 1 2010 ; 
I Watson et all 120091: IContursi et all 120071: iPeeters et al.l 12005 : 
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iPovich et alj|2007b . However, this phenomenon does not ap- 
pear to be systematic, and its origin is still poorly understood. 

In this paper, we present a systematic investigation of the 
properties of dust associated with the ionized gas i n the LMC . 
The LMC, located at a distance of ~ 50 kpc dFeaslI 119991: 
Ivan Leeuwen et all 12007). has an advantageous almost face- 
on viewing angle. At -32° latitude, the contamination along 
the LOS is expected to be sm all compared t o the plane of 
the Galaxy. For this galaxy, iBernard et al.l (120081) carried 
out a cross-correlatio n study of the combined Spitzer SAGE 
dMeixner et al.l 12006) and IRAS data, with HI and CO data. 
That work allowed to characterize (in terms of emissivities, 
abundances, etc.) the properties of dust mixed with the atomic 
and molecular medium, as well as to point out t he existence of 
an exc ess of emission at 70 fim. Interestingly, IBernard et al.l 
(2008) showed that this excess appears to be spatially cor- 
related with both the atomic and the diffuse (I# a < 200 R) 
ionized gas. 

The prese nt study build s over and expands the analysis 
described in Bernaix iet al.l (12008b . by including the ionized 
medium in the cross-correlation and by encompassing differ- 
ent physical regimes of this phase of the gas, from the most 
tenuous environments (diffuse medium), to the densest ones 
(HII regions). 

The paper is organized as follows. Section [2] summarizes 
the data we used. In Section[3] we describe the decomposition 
method as well as the dust emission modeling. We present our 
main results in Section [4] Concluding remarks are provided 
in Section[5] 

2. OBSERVATIONS 
2.1. IRdata 

Along the lines of Bern ard et al.1 (120081) . we use, as a tracer 
of dust emission from the NIR to the FIR, the Spitzer SAGE 
data. The data span the wavelength range 3.6 /mi to 160 fim, 
with an angular resolution from 2" to 40", respectively. 

These Spitzer data are combined with the 12 and 100 
fim IRIS data (Improved Reproces sing of the IRAS Survey, 
lMiville-Deschenes & Lagache 2005) characterized by an an- 
gular resolution close to 4'. The inclusion of the IRAC 8 fim 
and IRIS 12 fim data allows us to take into account the con- 
tribution from the most prominent PAH features, i.e. at 7.7 
and 1 1.3 fim, while with MIPS 24 and 70 fim the bulk of the 
emission from VSGs is incorporated. In addition, the 100 yum 
IRIS data point, combined with the MIPS 160 fim measure- 
ment, provides the necessary leverage for the derivation of the 
equilibrium temperature of the BGs. 

2.2. Gas tracers 
2.2.1. HI data 

As a tracer of atomic gas, we use the lKim et alJ d2003h 21- 
cm map, which is a combination of interferometric data ob- 
tained with the Australia Telescope Compact Array (ATCA), 
at 1' angular res olution, and the Parkes ante nna, at 15.3' an- 
gular resolution (iStavelev-Smith et al.l l2003). Assuming the 
gas is optically thin, the column density can be derived using 
the relation: 

N m = X HI W HI (1) 

where Whi is the HI integrated intensity map and the conver- 
sion factor Xh i is taken equal to 1.82 X 10 18 H/cm 2 /(Kkm/s) 
(Spitz er Hl978l) . The integrated map was built by summing 
up all the velocity channels from 190 to 386 km s . For the 



purpose of estimating and subtracting the Galactic foreground 
contribution to each IR wavelength , we use a Galactic HI col - 
umn density map constructed bv lStavelev-Smith et al.l 12 003). 
More details on this procedure can be found in Bernard et al. 
(120081 

2.2.2. CO data 

The 12 CO (J= 1-0) survey obtained with the 4-m radio 
NANTEN telescope of Nagoya University at Las Campanas 
Observatory, Chile (see Fuk ui et al.1 120081) is adopted as a 
tracer of the molecular gas. In particular, we make use of the 
second NANTEN survey which, characterized by an angular 
resolution of 2.60, evidenced th e existence in that galaxy of 
272 individual molecular clouds (Fukui et al. 20QH). 

We compute the molecular column densities using the rela- 
tion: 

N H2 = XcoWco (2) 

where Wco is the CO integrated intensity map, generated by 
integrating the original CO cube over the 200 < V] sr < 320 
km s _1 velocity range. Fuku i et al.l (120081) derived the Xco 
conversion factor assuming that each cloud is at virial equilib- 
rium. They found an average value for all the clouds of Xco = 
(7 + 2) x 10 20 H 2 /cm 2 /(K km s" 1 ). From cloud to cloud, val- 
ues vary from 1.6xl0 20 to 4.2xl0 21 H>/cm 2 /( Kkm/s). Their 
averag e value is larger than that estimated by Hughe s et al.l 
(2010) of 4.7 x 10 20 H 2 /cm 2 /(Kkm/s), using the Mopra tele- 
scope data, at a resolution of 33". This difference is due to 
the fact that low CO brightness clouds have been excluded 
from the MOPRA survey, and the Xco value has been de- 
duced from a smaller number of clouds compared to NAN- 
TEN observations. In this work, in order to b e consistent with 
our previous studies (e.g. Paradis et al. 2 01 ll) . we use the Xco 
values determined for each individual cloud by Fuk ui et al.1 
(2008). For clouds for which the virial mass could not be de- 
rived, we adopt the LMC average value of X C o=(7 + 2) X 10 20 
H2/cm 2 /(K km s~'). For a co mprehensive t abulation of Xco 
values, we refer the reader to Par adis et al.l (1201 ll) (see their 
Table 5). 

2.2.3. Ha and radio data 

The warm ionized gas, characterized by an electron tem- 
perature of the order of T e - 10 4 K, emits both recombination 
lines and radio free-free continuum. Among the observed re- 
combination lines, the Ha line is the brightest and is due to a 
transition of the Balmer line series at 6562.81 A. As for free- 
free emission, this consists of free electrons being scattered 
and decelerated by their encounters with ions. Ha and radio 
templates are available for the LMC and both tracers are used 
in our analysis to assure an independent cross-check of the re- 
sults. For the Ha , we use the Souther n H-Alpha Sky Survey 
Atlas (SHASSA. iGaustad et all 120011) data set which covers 
the southern hemisphere (6 < 15°), at an angular resolution 
of 0.8', and with a sensitivity of 2 Rayleigh/pixel (1 R=10 6 
photons cirT 2 s~'). Following Lag ache et al.l d 19991) . if we as- 
sume that the electron density n e is constant along each LOS, 
the H + column density can be derived using the relation: 

J!l_ = 1.37xl0 18 ^^ (3) 
Hem 2 R cm i 

where 1 R=2.25 pc c m -6 for T e =8000 K (see, for instance 
iDickinson et al.l 120031) . For n e , we adopt different values for 

5 about 40 pc at the distance of the LMC 
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diffe rent regimes of the ionized gas, as described in Section 
ED 

To trace the radio free-free emission, we use observations 
obtained with the P arkes 64-m telescope at 4.75 GHz (see 
Filipo vic et al.ll 19951 for a complete description of the Parkes 
survey). At that frequency the beam size is 4.9', but after cor- 
recting for scanning effects, the effective beam size results to 
be equal to 5.6'. 

In order to compare the radio and Ha maps and compute gas 
column densities for the radio map, we tra nsform the radio 
emissi on to Ha units, using Equation 1 1 in Dic kinson et alJ 
(2003): 



It 

Ilia 



8.396 x 10 3 a x Vgh' 7 ^ 667 ' 10 029/r4 (l + 0.08) 



(4) 



where Tb (brightness temperature) is in /jK and Iho. in 
Rayleigh. T4 is T e in units of 10 4 K, while the factor a de- 
pends o n the frequency and the electron temperature (see Ta- 
ble 3 in iDickinson et al.l 120031) and vghz is the frequency in 
GHz. For T e = 8000 K (T 4 = 0.8), a is close to unity. To 
apply Equation [4] we first convert the original units of the 
Parkes map (mJy/beam) into Jy, and these into Tb using the 
expression: 

_ S AJSGHzC 1 

with S4.75GHZ the total flux density at 4.75 GHz, c the speed of 
light, v the frequency and k the Boltzmann constant. Once the 
Parkes map has been converted into Ih» units, the H + column 
density can be computed by applying Equation|3] The sensi- 
tivity of the Parkes telescope is 8 mJy/beam, also equivalent 
to 28 R, i.e. significantly lower than the sensitivity threshold 
of the SHASSA map. This implies, as we will also see in 
Section 3.2, that while the SHASSA map will be used as a 
tracer of the ionized gas for all types of environmental condi- 
tions (i.e. from the diffuse to the dense ionized medium), the 
Parkes map will only be used to probe the densest/brightest 
regions. 

For the pixels corresponding to bright HII regions, we have 
performed a comparison of the radio Parkes map, converted 
into Ha emission (In ff ), with the SHASSA map. The com- 
parison highlights that, for 88% of the pixels, the Ha-Parkes 
data are higher than the SHASSA one. The median value of 
the ratio between the two data sets is 1.6. To explain such a 
discrepancy, we invoke extinction as a source of attenuation 
in the SHASSA map. In this respect, we remind the reader 
that Ha measurements are affected by dust absorption, while 
radio data provide an uncontaminated observational window 
on the ionized gas. In the LMC, extinction is usually con- 
sidered rather low. For instance, Ima ra & Blitzl (120071) and 
iDobashi et al.l ((2008) quote an average A(V) of 0.3 mag. In 
our analysis, the factor 1 .6 of discrepancy would induce an ex- 
tinction value in the SHASSA map of 0.5 1 mag in brig ht HII 
regions. Following |0'Donnelll(fT994l) . lDickinson et al.fi 2003) 
estimate the absorption in Ha as A(Ha) = 0.81A(V). This re- 
sult is obtained using a n extinction cu r ve for the Milky Way. 
However, as shown by Gordon et al. (2003), the extinction 
curve in the LMC is very similar to the Galactic one. There- 
fore, the estimate above also applies here. In that case, the cor- 
responding extinction of LMC bright HII regions in the visi- 
ble is A(V)=0.63 mag. This value is slightly higher than the 
ext inction estima t ed usi ng the LMC A(V) map constructed 
bv IDobashi etail d2008h . and weighted by the Ha emission 
Q^A(y)I Ha l YjiHa — 0.57 mag) over the same HII regions. 



It is unlikely that bright HII regions are all confined behind 
dust clouds, which may indicate that the extinction could have 
been underestimated when constructing the LMC A(V) map 
from star counts, as already suggested in Paradis et al. (201 1). 
In the following, we therefore expect differences of about a 
factor 1.6 when c ompar ing results using SHASSA and Parkes 
data (see Section l4TT3b . 



All the data used for the analysis (IR, HI, CO, Ha) have 
been degraded to the lowest angular resolution of the Parkes 
4.75-GHz map, using a Gaussian kernel with a fwhm of 5.6'. 
After being convolved, the maps have been reprojected onto 
the same grid, with a 2' pixel size. Finally, poi nt source sub- 
tractio n has been applied to all the IR maps (see lBernard et al J 
2008, for the description of the map processing). 

3. DECOMPOSITION OF THE IR EMISSION 
3.1. Regimes of the ionized gas 

We consider three regimes of the ionized gas, i.e. from dif- 
fuse H + to bright HII regions, and for each regime we in- 
vestigate the dust properties independently, due to potentially 
intrinsic variations in dust composition, evolution and physics 
of the grains. To define each regime, we set three thresholds 
in Ha emission intensity using the SHASSA map, convolved 
to 5.6' resol ution, and the catalog of HI I regions in the LMC 
compiled bv lKennicutt & Hod ge ( 1986) (hereafter KH86). 

We start by analyzing the pixel intensity distribution in the 
SHASSA map (see Figure[T]i. The mean (IhoO of the distribu- 
tion is 1 .4 R, with a dispersion (cr) of 1 . 1 R and a minimum 
of 0.4 R, which is the sensitivity of the SHASSA data for the 
5.6' angular resolution and the 2' pixel size of the map. About 
99% of the pixels lies in the range 0.4 R < Ih» < 6.7 R, cor- 
responding to l\i a + 5<t. We notice that none of the sources in 
the KH86 catalog are located in regions of the SHASSA map 
where the Ha intensity is within this range (see Figure|2]l, and 
therefore, we adopt it to define the diffuse ionized gas regime. 

For l\ia > 6.7 R, we look, in the SHASSA map, at the inten- 
sity distribution of the pixels which correspond to the location 
of the HII regions in the KH86 list. To select the pixels, we 
use, for each cataloged source, its coordinates and measured 
angular diameter. The Ha intensity distribution obtained by 
this process is shown in FigureQ~](insert panel). This is charac- 
terized by a mean value of 26.3 R, with a dispersion of 17.4 R. 
In this case, 99% of the pixels lies in the range 6.7 R < Ih» < 
1 13.3 R, and such a range is used to define what, in this con- 
text, we denote as typical HII regions. 

For Ih» > 113.3 R, we have a few very bright HII regions 
which are listed in the KH86 catalog. We set an upper limit 
of 1000 R to exclude extremely bright sources, such as 30 
Doradus, which - due to their complexity - would require a 
dedicated analysis. The range 113.3 R < In a < 1000 R hence 
identifies the third regime. 

In summary, we have three distinct regimes/cases which we 
will be tackled by our analysis: 

• case 1: 0.4 R < lu a < 6.7 R (diffuse ionized gas); 

• case 2: 6.7 R < Ih» < 113.3 R (typical HII regions); 



case 3: 113.3 R < I Ha 
gions). 



1000 R (very bright HII re- 



For each regime, we need to adopt a characteristic n e in or- 
der to compute the H + column density from Equation [3] For 



TABLE 1 

Results of the IR emission decomposition with respect to different gas phases, in unit of MJy/sr for 10 20 H/cm 2 , and their 1<t uncertainties. The different regimes of the ionized gas described in Section |3~T1 

ARE DENOTED - IN THE FIRST COLUMN - WITH 1 (DIFFUSE IONIZED GAS), 2 (TYPICAL HII REGIONS) OR 3 (BRIGHT HII REGIONS; 3A-USING THE Hff DATA AND 3B-USING THE RADIO DATA AT 4.75 GHz). 
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-3 
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the diffuse component, no estimate is available from the liter- 
ature. We have made an attempt to derive this value directly 
from the SHASSA map. For this purpose, we first compute, 
for each pixel in the range 0.4 R < < 6.7 R, the emission 
measu re (EM), by applying Equation |4] in Dic kinson et al.l 
(2003) and assuming T e = 8000 K. To derive the rms elec- 
tron density from the average EM so computed (EM = 5.3 + 
0.1 pc cm -6 ), we need to assume a linear scale of the emit- 
ting ionized medium. Our 5.6'-SHASSA map has a pixel size 
of 2'. At a distance of 50 kpc, this corresponds to a linear 
size of ~ 30 pc. However, such a value is very likely a se- 
vere underestimate of the actual vertical extent of the diffuse 
ionized gas. In fact, KH86 have analyzed the distribution in 
linear diameters for their cataloged HII regions, and obtained 
that the largest sources in the sample have a size of about 400 
pc. In the hypothesis that the diffuse gas extends far beyond 
the layer of discrete HII regions, if we denote with de+ the 
vertical extent of H + and take d H + = 500 pc, we obtain an av- 
erage electron density n e , of 0.10 + 0.08 cm 4 . This value is 
close to what has been derived for our Galaxy, i.e. n e ~ 0.03 
- 0.08 cm -3 , for a n estimated scale h eight ranging between 
1000 and 1800 pc (Haffn er et al.ll2009l) . If, instead of 500 pc, 
we assume de+ = 1000 pc, we retrieve n e , = 0.07 + 0.06 cm" 3 , 
matching the Galactic value. Despite such a result, given the 
uncertainty in de+ for the LMC, we choose to adopt the av- 
erage value derived from Galactic measurements. Therefore, 
for case 1, we take n e , = 0.055 cm 4 . 

For cases 2 and 3, we compute EM, for each HI I region in 
the KH 86 catalog by applying Equation 4 in Dic kinson et alJ 
(2003). For each source, we use the quoted Ho- flux and we 
set T e = 8000 K. The measured Ha flux is converted into 
Rayleighs, and the intensity cuts at 113.3 R and 1000 R for, 
respectively, case 2 and 3 are applied. The averaged emission 
measures for the two cases (EM2 = 183.8 + 17.8 pc cm 4 , 
EM3 = 1322.4 + 167.7 pc cm 4 ) are then used to derive an 
rms electron density. Following this procedure, we obtain: 
n e2 = 1.52 ± 0.44 cm 4 and n e , = 3.98 + 2.05 cm 4 , i.e. 20 
and 80 times larger than the diffuse component. 

3.2. Method 

In the LMC, iBernard et al.l (120081) evidenced the existence 
of a large FIR excess with respect to the atomic, molecular 
and ionized gas phases, which they attributed to either molec- 
ular gas with no associated CO emission or to optically thick 
HI. They showed that t his co mponent is spatially correlated 
with HI. Paradis et aTj (1201 ll) performed a statistical analy- 
sis of LMC molecular clouds and showed that if the FIR ex- 
cess was not taken into account in the IR correlations, the 
dust abundance in the atomic phase would be systematically 
overestimated by a f actor 1.7 (see their Section 3). Following 
Paradi s~et al.l (1201 11) . we decompose the IR emission at each 
wavelength and for each regime, as a function of the atomic, 
molecular and ionized gas column density: 

1(A) = a(Nf + N§) + bN^° + cN% + + d (6) 

with the column density in the FIR excess component. In 
the following, we refer to the sum of the atomic and FIR ex- 
cess components as the atomic phase. Despite this working 
hypothesis, it is important to keep in mind that the FIR excess 
could be molecular in na ture, instead of atomic, and associ- 
ated with warm H2 gas (Planck collaboration 201 le). The 
parameters a, b, c, denoting the emissivities associated with 
each phase of the gas at a given wavelength, and the constant 



d are determined using regression fits. To remove any possible 
offsets in the data, we subtract a background in all the maps. 
The background is computed as the median over a common 
area, corresponding to the faintest 10% of the HI data. This 
step also ensures to have a null IR emission for a null column 
density. The Ho- data are used for all regimes of ionized gas 
described in the previous section, while, due to the limited 
sensitivity of the 4.75-GHz map, we can only use the radio 
data for case 3. As a consequence, a direct cross-check of the 
results obtained by adopting both tracers will only be possible 
for this case. 

Although case 3 corresponds to bright HII regions for 
which, presumably, most of the radio emission is due to ther- 
mal free-free, yet we might have a residual non-thermal con- 
tamination associated with synchrotron radiation. To investi- 
gate the synchrotron contamination, we complement the 4.75- 
GHz Parkes data with the 1.4-GHz data, obtained with the 
same telescope. Both maps are smoothed to the common res- 
olution of 16.6', ch aracterizing the 1 .4-GHz data. Following 
iNiklas et al.l d 1997b . we estimate the fraction of emission con- 
tributed by free-free as: 

^=/.(^r + (i-A)(^p 

S 1.4 \VlA ) \ v 1.4/ 

where f t h is the thermal fraction at 4.75 GHz, S4.75 and S1.4 are 
the total flux densities at 4.75 and 1 .4 GHz, and as and a sync 
are the free-free and sync hrotron spectral ind ices. Following 
iPaladini et all ( 120051) and lHughes et all d2006l) . we take a ff = 
0. 1 . The spectral index of synchrotron radiation is known to 
present large spatial variations. In our Galaxy, it ranges be- 
tween 0.6 in star-forming complexes to 1 .2 in the interarm re- 
gions. If, in Equation|7] we adopt a fixed value of a sync =0.6, 
we derive an average thermal fraction in the LMC of 0.85. 
This value is probably an underestimate due to the fact that, 
while QTff = 0.1 is a good approximation for frequencies less 
than or equal to a few GHz and T e = 8000 K, the actual free- 
free spectral index is a weak function of both frequency and 
electron temperature (Be nnett et al.ll 1992b . Moreover, HII re- 
gions span a wide range of T e values. If we compare our re- 
sult with simila r attempts f o und in the literature, we see that, 
for instance, iMeinert et al.l (119931) tried to separate the ther- 
mal and non thermal component in the radio emission of the 
LMC, using Ha and FIR data, at a resolution of 51'. At 4.75 
GHz, they obtained a thermal fraction of 0.59. However their 
result applies to both the diffuse medium and discrete HII re- 
gions, while we only consider HII regions. 
In order to take the synchrotron contamination into account, 
we multiply the 4.75-GHz radio map by the derived thermal 
fraction equal to 0.85. 

Results of the correlations for each phase and regime of the 
gas are presented in Figure [3]and in Table Q] 

3.3. Modeling 
3.3.1. Single radiation field along the LOS 

From Equation|6]we obtain the emissivities associated with 
each phase of the gas for each wavelength A in the 3.6 fim 
- 160//m range, and with these values we build the corre- 
sponding spectral energy distributions (SEDs). The SE Ds are 
modeled with DustEM dCompiegne et all 120081 120 1 01) The 
DustEM code is an updated version of the lDesert et al.f(ll990l) 
model. The code allows to constrain mass abundances (Ybg, 
YvsGj and YpahX relative to hydrogen, of the three dust popu- 
lations. The intensity of a NIR continuum, with a temperature 
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TABLE 2 

Results of the modeling obtained by assuming a single radiation field along the LOS. Three different regimes of the ionized gas are considered: diffuse (1), 

TYPICAL HII REGIONS (2) AND BRIGHT HII REGIONS (3a CORRESPONDS TO Hff DATA AND 3b TO RADIO DATA AT 4.75 GHz). TH E FIRST COLUMN DENOTES THE ADOPTED RADIATION 

field spectrum: the Mathis RF or the 4-Myr cluster Galev RF (see Section |3.3.U . 



RF 


Cases 


Y to , (DIG) 


NIR cont 


Xrf 


Ypah/Ybg 


Yvsg/Ybg 






10- 3 


1(T 4 






10- 2 


Mathis 


1-HI 


1.82±0.04 


1.98±0.22 


1.51 ±0.03 


4.29±0.20 


5.37±0.28 


Mathis 


1-CO 


3.19±0.87 


4.90±1.36 


0.61±0.16 


9.87±5.19 


6.65±3.78 


Mathis 


1-Ha 


0.50±0.07 


36.71 + 1.83 


6.46+0.90 


1.50±0.54 


9.47±2.80 


Mathis 


2-HI 


1.62±0.03 


1.02±0.15 


1.86±0.03 


4.00±0.14 


7.08±0.26 


Mathis 


2-CO 


2.77±0.17 


14.44±0.69 


1.42±0.08 


8.94±1.07 


10.55+1.35 


4-Myr cluster 


2-Ha 


3.16±0.15 


91.69±3.83 


0.39±0.02 


0.88±0.09 


6.88±0.66 


Mathis 


3a-HI 


1.29±0.08 


3.06±0.33 


3.93±0.24 


3.27±0.41 


13.81+1.90 


Mathis 


3a-CO 


3.61±0.35 


16.47±1.07 


2.79±0.26 


6.38±1.20 


9.13±2.10 


4-Myr cluster 


3a-Ha 


12.67±0.94 


157.06±4.41 


0.21± 0.02 


0.77±0.12 


17.35+2.61 


Mathis 


3b-HI 


1.10±0.10 


0.78±0.37 


3.66±0.30 


3.54±0.59 


10.28±2.08 


Mathis 


3b-CO 


3.34±0.36 


13.77±1.12 


2.65±0.28 


6.83±1.45 


6.26±1.95 


4-Myr cluster 


3b-Radio (5GHz) 


7.33±0.51 


87.89±2.65 


0.22±0.01 


().68±0.10 


15.48±2.20 



TABLE 3 

Results of the modeling obtained by assumtng a composite radiation field along the LOS. Three different regimes of the ionized gas are considered: diffuse 

(1), TYPICAL HII REGIONS (2) AND BRIGHT HII REGIONS (3a CORRESPONDS TO Ha DATA AND 3b TO RADIO DATA AT 4.75 GHz). THE FI RST CO LUMN DENOTES THE ADOPTED 

composite radiation field spectrum: the Mathis RF or the 4-Myr cluster Galev RF (see Section |3.3.2> 



RF 


Cases 


Y, ot (D/G) 


a 


Ypah/Ybg 


Yvsg/Ybg 






io- 3 




IO- 2 


IO- 2 


Z, Zj Mathis 


1-HI 


3.39±0.02 


2.47 


4.57±0.06 


2. 11 ±0.08 


£i Zj Mathis 


1-CO 


2.35±0.13 


2.50 


8.68±0.78 


4.17±0.91 


£i Zj Mathis 


1-Ha 


3.06±0.12 


2.29 


2.74±0.28 


0.43±0.50 


Zi Zj Mathis + 4-Myr cluster 


2-HI 


6.23±0.03 


2.50 


2.42±0.03 


1.56±0.05 


Zi Zj Mathis + 4-Myr cluster 


2-CO 


7.83±0.11 


2.50 


4.91±0.12 


3.93±0.22 


Zi Zj 4-Myr cluster 


2-Hff 


1.81±0.02 


2.50 


0.88±0.03 


1.02±0.09 


Zi Zj Mathis + 4-Myr cluster 


3a-HI 


8.38±0.16 


2.32 


1.66±0.07 


0.27±0.26 


Zi Zj Mathis + 4-Myr cluster 


3a-CO 


19.38±0.58 


2.50 


3.38±0.17 


3.81+0.55 


Zi Zj 4-Myr cluster 


3a-Hff 


3.53±0.08 


2.50 


0.76±0.04 


8.39±0.39 


Zi Zj Mathis + 4-Myr cluster 


3b-HI 


7.17±0.19 


2.38 


1.89±0.10 


0.72±0.41 


Zi Zj Mathis + 4-Myr cluster 


3b-CO 


17.01+0.59 


2.50 


3.61±0.20 


2.08±0.64 


Zi Zj 4-Myr cluster 


3b-Radio (5GHz) 


2.16±0.04 


2.50 


0.71±0.04 


7.00±0.36 



of 1000 K, is also a free parameter of the model. Indeed, the 
inclusion of this component is required, as discussed in sev- 
eral studies (e.g. iLu et al.ll2003l) . and allows us to reproduce 
observational data at 3.6 and 4.5 /an. Two different radiation 
fields (RFs) input models are used in this work. The strength 
of the field (Xrjj) is a free parameter in both cases. How- 
ever, for all 3 cases (diffuse medium, typical and very bright 
HII regions), we impose, for dust emission associated with 
the atomic and molecular phases, as well as for dust emis- 
sion associated with the diffuse ionized gas in case 1, the 
same interstellar R F hardness as for the solar neighborhood 
(Mathi s et aT]|1983l) . In the following we will refer to this RF 
as Mathis RF. We note that, for dust in the atomic and molecu- 
lar gas phases, we assume that the dominant component of the 
RF is always the one related to the interstellar solar neighbor- 
hood. However for cases 2 and 3, dust in the ionized phase is 
associated with typical and bright HII regions, and therefore 
it is essentially heated by local hot stars. In these regions, the 
solar neighborh ood interstellar RF is no t a good approxima- 
tion. Following Ka wamura et alJ (|2009), the youngest stellar 



objects are less than 10 Myr. Using the GALEV code0, which 
is an implementation of evolutionary synthesis models, we 
have generated an UV/visible spectrum of young stellar clus- 
ters of 4 Myr. In the following, we will refer to this RF as the 
4-Myr cluster RF. The summary of the RFs used for each case 
and phase is described in Table [2] and the RF spectral shapes 
are shown in Figure |4] 

3.3.2. Composite radiation fields along the LOS 

Assuming a single RF component along the LOS is a sim- 
plistic hypothesis. Therefore we tested our results using in- 
stead a combination of several RF intensities and hardnesses 
along the LOS. For case 1, we apply the "loca l-SED com- 
bination" model proposed by iDale et aD (120011) for normal 
star forming galaxies. This model assumes a power-law dis- 
tribution of dust mass subject to a given heating intensity 
dM d (XRF): 

dM d (X RF ) oc X R a F dX RF , (8) 



6 available at http://www.galev.org 
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Fig. I. — Pixel intensity distribution of the SHASSA Ha map. Values on the x-axis have been truncated to 30 R to evidence the bulk of the distribution. Overlaid 
are: (solid line) the mean (xo) and (dashed line) 5<r values. Insert panel: pixel intensity distribution for the KH86 HII regions in the SHASSA map. The mean 
value and 5<r limits are also shown. 
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Fig. 2. — SHASSA Her map (a) and spatial locations of pixels (in grey) selected for each case: (b)- case 1 (diffuse ionized gas), (c)- case 2 (typical HII regions), 
and (d)- case 3 (bright HII regions). The overlaid symbols show the HII regions (Kennicutt & Hodge 1986). The range of the maps is [-8 ; 100] R. 



where a controls the relative co ntribution of the various 
strengths of the field to the SEDs. iDale et all (12001b consid- 
ered RF strengths and a values in the range 0.3 < Xrf < 10 5 
and 1 < a < 2.5, respectively. A value of a close to 1 in- 
dicates an active star-forming region, whereas a around 2.5 
corresponds to the diffuse medium. This parameter is linked 
to the Igo/Iioo ratio, and therefore gives an indirect measure of 
dust temperature. 

Here is how we proceed in our analysis. We first compute 
emission spectra using DustEM, for different Xrf values, and 
the Mathis RF (noted I™ d (X RF , RF )). Then we sum up these 
emission s pectra over the same Xrf range as proposed by 
IDale et al.ld2001l) . according to: 



Tii H jl'v" (Xrfj, RFq) X X R 
2/ HjXrfj 



(9) 



For the atomic and molecular phases in cases 2 and 3, we 
check the impact of adopting a combination of the Mathis RF 



and the simulated 4-Myr cluster Galev RF. Equation [9] be- 



jtOl _ 



2 2i Yij I'v" (XrfJ, RF Q ) X X RF j 
\f Hi Z j F"" d (X RFi , RF Galev) x 



(10) 



RF.i 



2/ YjjX RFi 

The coefficient / allows to normalize the emission spectra 
computed with each RF hardness. Indeed, the same Xrf will 
not give the same absolute emission level depending on the 
RF used. To determine /, we adjust the emission spectra using 
the two RFs (Mathis RF and 4-Myr cluster RF) and look for 
the coefficient allowing the same brightness at the maximum 
emission. This coefficient / is re-derived for each Xrf. For 
the ionized phase in cases 2 and 3, we apply Equation [9] but 
replacing the Mathis RF by the 4-Myr cluster RF. 
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The DustEM code does not allow to perform a model fit us- 
ing a combination of RFs. Therefore, we have to use a differ- 
ent method to analyz e the effect of a composite RF. Following 
iParadis et al.l (l2009al) . in order to derive mass abundances for 
each dust grain species, we pre-compute the IR brightnesses 
associated w ith each gas phase, using the Galactic abundances 
obtained by iDesert et al.l dl990l) . for different a values. We 
then compare the observed brightnesses between 8 and 160 
/mi with the predicted values, and infer the dust mass abun- 
dances minimizing the \ 2 f° r eac h value of a. We select the 
solution with the lowest^ 2 . The additional NIR continuum is 
not taken into account to compute the IR brightnesses, since 
we do not know the origin of this component, nor how to treat 
it in the case of a combination of RFs. For this reason, we do 
not consider IR data below 8 pm. The summary of the RF 
combinations used for each case and phase is given in Table 
® 

4. RESULTS AND DISCUSSION 

The dust emission spectra are shown in Figure [3] Figure 
[5] illustrates the SED modeling obtained using DustEM, as- 
suming a single radiation field component. In the following 
we analyze the results for the diffuse ionized gas, typical and 
bright HII regions. We also discuss their consistency using the 
hypothesis of a RFs combination. The best-fit parameters for 
each phase of the gas (molecular, atomic and ionized), and 
each regime of the ionized gas, as well as for different RFs 
assumptions, are provided in Tables [2]and[3] 

4.1. Single radiation field component 
4.1.1. Case 1: diffuse ionized gas 

Figures [3] (panel A) and [5] indicate that, when compared 
to the atomic and molecular phases, dust emission associ- 
ated with the ionized gas phase is a relevant contribution 
to the total emission. However, we emphasize that the 
absolute amplitude of the spectrum for the ionized phase is 
proportional to n e . As a consequence, if n e has been under- 
estimated/overestimated, the absolute level of the spectrum 
has also been underestimated/overestimated, although the 
spectral shape will remain unchanged. Here, as discussed 
in Section 3.1, we take n e equal to 0.055 cm -3 , which is 
the average observe d value in our Galaxy, as reported by 
lHaffneretal.l(l2009l 

In this case, the adopted RF hardness is the same for all gas 
phases, as discussed in Section[33J] Therefore, it is possible 
to perform a direct comparison among the temperature esti- 
mates resulting from the best-fit RF intensities. We remind 
the reader that, if we assume a greybody dust emission model 
with spectral index then the RF intensity and the equilib- 
rium temperature of the BGs are related through: 



X 



RF 



RF 



T LMC \ 
l _d_ 

1 d 



4+/J 



(ID 



where X° F 
(X° F - 1), 
hood is 
with p = 2. 
K and 23.9 
respectively, 
the ionized 
molecular g 



is the RF scaling factor for the Mathis RF 
and the dust temperature in the solar neighbor- 
taken equal to 17.5 K dBoulanger et alj IT996). 
We find dust temperatures equal to 18.7 K, 16.1 
K for the atomic, molecular and ionized phases, 
These results confirm that dust associated with 
phase is warmer than dust in the atomic and 
as phases, although our derived temperature for 



the diffuse H + is sli ghtly lower than the value obtained by 
iPaladini et all d2007l) . of 26.6 + 0. 1 K. As for dust emissivity 
for the ionized phase, we obtain an e stimate of 2.3x 10~ 26 
cm 2 /H at 160 pm. For comparison, Laga che et al.l (120001) 
find, at high latitude, an emissivity of l.lx 10~ 25 cm 2 /H at 
250 pm, which corresponds to an emissivity of 2.7x 10~ 25 
cm 2 /H at 160 pm ass uming a spectral index equal to 2, while 
Paladi ni et al.l (|2007), for the Galactic plane, report a value 
of 35.1 MJy/sr for 10 2l) H/cm 2 at 140 pm, corresponding to 
8.7 x 10~ 25 cm 2 /H at 160 pm. Therefore, even if in terms of 
temperature, dust in the ionized phase of the LMC is roughly 
comparable to its Galactic plane counterpart, it appears to 
be 40 times less emissive in the LMC with respect to the 
Galactic plane, and 10 times less emissive with respect to the 
high Galactic latitudes. The low value in the LMC cannot 
be only due to the low m etallicity of the LMC (Z ^ 0.3 
- 0.5 Z Q [Westerlund] 119971) with respect to our Galaxy, as 
the discrepancy in emissivity is significantly larger than the 
difference in metallicity. For the Galactic plane, at least part 
of the difference could be induced by confusion along the 
LOS. However, we think that different dust properties in the 
ionized gas phase of the LMC compared to our Galaxy (high 
latitudes and Galactic plane) are likely at the origin of the 
observed discrepancy. 

From Table|2] the total amount of dust in the ionized phase, 
similar to the dust-over-gas ratio (D/G), appears to be 3.6 and 
6.4 times lower than in the atomic and molecular phases. In 
particular, the PAH relative abundance with respect to BGs 
reveals a decrement by a factor 2.9 and 6.6 when compared 
to the atomic and molecular phases, respectively. The high 
value of PAH relative abundance in the molecular phase is 
connected to the low Xrf which, in turn, is a consequence of 
the shielding effect from the RF induced by the higher den- 
sity characterizing this phase of the gas. For the VSG relative 
abundance in the ionized phase, we do not find the same de- 
pletion behaviour as for PAHs. The VSG relative abundance 
is in fact roughly the same in the three phases. We interpret 
this effect as VSGs being resiliant to the RF. 

Remarkably, the spectrum of dust associated with the ion- 
ized phase highlights more emission in the NIR domain, by a 
factor of 7.5 and 18.5, compared to the mole cular and atomic 
phases. The origin of this effect is not clear. iLu et al.l (120031) 
argue that the NIR continuum intensity is quite linear with 
the strength of the aromatic features in emission. They sug- 
gest that both this continu um and the aromati c features origi- 
nate from similar carriers. FlagexetalJ(|2006), using Spitzer 
data for the Galactic diffuse interstellar medium, require the 
presence of a continuum to explain the observations, but they 
conclude that the continuum carriers might not be the PAHs. 
What we observe in the ionized phase of the LMC is a de- 
crease of the PAH relative abundance and an increase of the 
NIR continuum, whereas the low value of the NIR contin- 
uum intensity in the molecular phase is associated with a high 
val ue of PAH relative abundance. Therefore our results favour 
the lFlagev et all d2006V s interpretation. 

4.1.2. Case 2: typical HII regions 

In this case, the modeling seems to reveal (Figure [3] panel 
B) that emission associated with the ionized gas phase is not 
only a major contributor to the total emission, but is even 
dominant with respect to the other components, i.e. atomic 
and molecular. Nevertheless, as in the previous section, we 
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Fig. 3. — Spectral energy distribution associated with the atomic (solid line), molecular (dashed line) and ionized (dashed-dotted line) phase of the gas, for case 
1 (diffuse ionized gas, panel A), case 2 (typical HII regions, panel B) and case 3 (bright HII regions). For case 3, we show the results of the analysis obtained by 
using both Ha (case 3a, panel C) and 4.75-GHz data (case 3b, panel D). Panels E, F and G: spectral energy distributions of all cases (case 1 in thin black line, 
case 2 in thin grey line, case 3a in thick black line and case 3b in thick grey line), plotted at once for each gas component (same notation as in panels A, B, C and 
D). 



need to remind the reader that this result quite strongly de- 
pends on the electron density. At the same time, we do not 
think that the adopted value can be off by a large amount, 
given that it has been derived from the HII regions in the 
HK86 catalog. In fact, we should be wrong by a factor 5-6 
in order for the ionized gas spectrum to match the atomic and 
molecular spectra at 160 /mi. 

For the ionized phase associated with typical HII regions, 
we adopt - as discussed in Section 13.3.11 - a bluer RF input 
model. Indeed, the Mathis RF is not representative in this 
case, as the ionized gas is mostly heated by radiation from 
nearby stars. Since we assume different RFs for different gas 
phases, we cannot directly compare the temperature inferred 
from the model. However, Figure |5] evidences a shift of the 
BG emission peak towards shorter wavelengths in the ionized 
gas phase compared to the other phases. This result would 
suggest an increase of dust temperature when going from the 
atomic or molecular phase to the ionized one. The hardness of 



the RF modeled with Galev is proportional to A . We have 
checked our results against different slopes of the input RF, 
for example A and /T 15 (see Figure [4j, and verified that 
the changes in the outputs are negligible. 

As for case 1 , the ionized phase appears to be characterized 
by less PAHs relative to BGs. In particular, the PAH relative 
abundance decreases by a factor 4.5 to 10, compared to the 
relative abundance in the atomic and molecular phases. This 
is likely a consequence of a higher RF intensity with respect 
to the other phases. 

Furthermore, in a similar fashion to what we find for the 
diffuse ionized medium, the molecular phase presents an en- 
hancement of the PAH relative abundance. For typical HII re- 
gions, the CO data trace the photo-dissociation region (PDR) 
surrounding the ionized gas, where the local high density con- 
tributes to shield the PAHs from the stellar RF. Interestingly, 
the atomic phase also shows a decrease in the PAH relative 
abundance compared to the molecular phase, even if the ef- 
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Fig. 4. — Radiation field templates used in this work: 



Mathis (solid black line), 4-Myr cluster Galev (solid grey line), 



removed) with a hardness proportional to A 03 and A 15 (dashed black and dashed grey lines). 



and 4-Myr duster Galev (gas lines 



feet is less important than for the ionized phase. The atomic 
phase is subject to an attenuated stellar RF, as well as to mod- 
erate radiation from the ISM. The combined effect of both 
RFs might cause the partial PAH destruction, or a decrease of 
the excitation state of the PAHs. In this latter scenario, the 
PAHs would be less emissive. Overall, the results obtained 
for typical HII regions confirm those derived for the diffuse 
H + . The VSG relative abundance is qualitatively the same in 
the three phases. 

As for case 1, we notice that the NIR continuum intensity 
is increased in the ionized phase with respect to the atomic 
phase, in this case by a factor 90, whereas the PAH relative 
abundance decreases. This result corroborates the hypothesis 
that the NIR continuum is not correlated with the PAH abun- 
dance. 

4. 1 .3. Case 3: bright HII regions 

In case 3, we are looking at dust emission in the vicinity 
of conspicuous massive star associations. As in the previous 
case, dust emission in the ionized phase appears to dominate 
the total emission (see Figure [3] panels C and D), although 
with the caveat discussed in Section 4.1.1 and 4.1.2. 

Results derived using Ha and radio data are in good agree- 
ment. The lower amplitude of the spectrum associated with 
the ionized phase when in Equation [6] we make use of the 
Parkes rather than the SHASSA data is likely partly a con- 
sequence of dust absorption in the SHASSA map, as already 
mentioned in Section |2.2.3l 

From Table [2] we notice that the RF intensity and the 
dust relative abundances are in good agreement for both 
the ionized gas tracers. The results show, for the ionized 
gas phase in particular, a low PAH relative abundance 
comparable to case 2, and an increase of the VSG relative 
abundance by a factor 2.2-2.5 (Yvsg/Ybg — 6.9 in case 
2 and Y V sg/Ybg - 15.5 - 17.4 in cases 3). Remarkably, 
these conclusions do not depend on the specific spectral 
shape of the RF modeled with Galev. The molecular phase 
appears again to be characterized by a high PAH relative 
abundance. We note that the PAH destruction in the ionized 
phase becomes more pronounced going from case 1 to 
case 3, probably due to the increase of the RF intensity 
from the diffuse medium to bright HII regions. This is 
consistent with previous findings repo rted in the lit e rature . 
For instance, using ISOCAM data, IContursi et al.l (2007) 



found a lack of aromatic features carriers in the core of the 
LMC HII region N4, and suggested that the high and hard 
RF might be respon sible for the PAH destruction. Likewise, 
iPovich et all (120071) . id entified PAH destructio n at the edge 
of the HII region M17. iBoulanger et alj dl988l) explained the 
decrease of the 12 /mi emission in the California Nebula by 
grain destruction in regions with high UV radiation. They 
estimated that around 80% of the 12 yum emitters could be 
destroyed whenever the energy density is larger than 50 times 
the energy density in the solar neighborhood. Despite this 
theoretical and observational evidence, we cannot state that 
the RF is systematically responsible for the depletion of 
PAHs in HII regions. Indeed, several other processes could 
be at work, like Coulomb destruction of PAH dications or 
de hydrogenation by dire ct photodissociation, as suggested 
by iBernard et al.l ([T994) based on the analysis of the 3.3 
pm emission and nearby continuum in the planetary nebula 
BD+30 3639. Moreover, iGiard et al.l (119941) studied the same 
dust feature in the Orion bar and M17, and propose that H + 
chemical sputtering could cause PAH destruction. 

In terms of temperature, as in case 2, we can deduce from 
Figure [5] the presence of warmer dust in the ionized gas as 
compared to the other phases, from the fact that the peak of 
emission for the BGs is located at shorter wavelengths. Fi- 
nally, as for case 1 and 2, the NIR continuum is higher in the 
ionized phase than in the other phases. 

4.2. Composite radiation fields 

In this section, we discuss the use of a combination of RFs, 
instead of a single RF, for fitting the SEDs corresponding to 
each phase of the gas (see TableO. 

For case 1, the values given in Table [3]confirm the decrease 
of the PAH relati ve abu ndance in the ionized phase, as evi- 
denced in Section l4.1.1l On the contrary, the enhancement of 
the PAHs in the molecular phase might suggest that these en- 
vironmental conditions are indeed particularly favourable for 
their survival. For the VSG relative abundance in the ionized 
phase, we find a significantly lower value with respect to what 
we obtain by adopting a single RF component (see Tables [2] 
and [3]). However, the VSG relative abundance is highly de- 
pendent on the quality of the fit at FIR wavelengths and, in the 
ionized phase, the model exhibits an important lack of emis- 
sion at 70 yt/rn, as well as a higher^ 2 , as compared to a single 
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Fig. 5. — Results of the fits obtained with the dust emission model (DustEM) assuming a single RF along the LOS, for each gas phase and for each regime of 
the ionized gas : diffuse (case 1), typical HII regions (case 2) and bright HII regions (case 3). In particular: case 3a corresponds to Ha data; case 3b to radio data. 
The squares show the fits of the model after applying the color-correction in each Spitzer and IRIS bands. Different curves denote the contributions from various 
grain populations: total (solid), BG (dot-dash), VSG (dot), PAH (dash), and the NIR continuum (long dash). 



RF modeling. This might indicate that the low VSG relative 
abundance has, for a composite RF, been underestimated. 

For case 2, the combination of the Mathis with the 4-Myr 
cluster RFs, both for the atomic and molecular phases, gen- 
erates a bluer RF than a single Mathis RF. Therefore, as ex- 
pected, we find a decrease in the PAH and VSG relative abun- 
dances. Indeed, PAHs and VSGs absorb most of their energy 
at short wavelengths, while BGs absorb energy over the whole 
RF spectrum. As a consequence, PAHs and VSGs absorb and 
re-emit more energy with a bluer RF, and to produce the same 
IR flux with a bluer RF than the Mathis RF, the model requires 
less PAHs and VSGs. For the ionized phase, the combination 
of 4-Myr cluster RFs gives the same PAH relative abundance 
of a single 4-Myr cluster RF, whereas a significant decrease of 
the VSG relative abundance is observed. However, the gen- 
eral trend is similar to what found for a single RF: as in case 1 , 
the highest PAH relative abundance is in the molecular phase, 



which also exhibits more VSGs with respect to BGs, while 
the atomic phase shows intermediate values of PAH and VSG 
relative abundances, compared to the other gas phases. 

In case 3, as for case 2, the PAH relative abundance in the 
atomic and molecular phases decrease by a factor ~ 2 com- 
pared to a single RF. In the ionized phase, we still have less 
PAHs than in the atomic phase, by a factor 2.2 and 2.7, for 
case 3a and 3b, respectively, and by a factor between 4.4 and 
5.1 compared to the molecular phase. Results are therefore in 
agreement with the modeling with a single RF. 

For the VSG relative abundance, we obtain, for the atomic 
phase, different and contradictory results depending on the 
assumption on the RF, therefore preventing a clear interpre- 
tation of the results. In the ionized phase, the modeling con- 
firms the enhancement of the relative VSG abundance already 
highlighted by a single RF. 

Overall, the results obtained assuming a composite RF sub- 
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stantially corroborate those derived with a single RF. 

4.3. Evolution of dust emission associated with different gas 
components 

In panels E, F and G of Figure [3] we present the results 
of the analysis in a slightly different fashion, i.e. by com- 
bining in individual plots, each corresponding to a given gas 
phase, the SEDs obtained for all cases (1,2 and 3). We can 
see that no matter which gas component we consider, the am- 
plitude of the spectrum always increases going from the dif- 
fuse medium to bright HII regions. However, for the atomic 
phase, the difference in amplitude between the different cases 
is less important compared to the other gas components, and 
a rather remarkable agreement among cases is found at 8 fim. 
As previously noted, the absolute level of dust emission in 
the ionized gas is directly proportional to n e , for which the 
adopted value changes according to the specific environment. 
Therefore, the discrepancy in absolute amplitude of the spec- 
tra for that phase is not well constrained. On the contrary, for 
the atomic and molecular phases, we are confident that the re- 
covered SED amplitudes are reliable, as these depend on the 
HI and CO conversion factor (Equation 1 and 2) which, re- 
spectively, is known to be constant for a given environment 
(HI), and has been accurately determined for the majority of 
the clouds (CO). 

We observe prominent 8 yum emission in the atomic phase, 
probably associated with the 7.7 p.m PAH feature. This emis- 
sion feature does not show in the ionized phase for the case 
of HII regions, likely due to PAH depletion. In comparison, 
the molecular phase is clearly characterized by strong 8 fim 
emission, as well as by significant emission at 5.8 and 12 fim. 

In the FIR domain, the atomic phase presents an excess 
of emission at 70 fim, as does the ionized ph ase in the dif- 
fuse m edium. This finding is in agreement with B ernard et al.l 
J2008h . 

5. CONCLUSIONS 

Combining IR data with tracers of HI, H2 and H + , we have 
evidenced, for the first time, dust emission associated with the 
ionized gas in the LMC. We have performed our analysis for 
different regimes of the ionized gas : diffuse gas (case 1), typ- 
ical HII regions (case 2) and compact/bright HII regions (case 
3). We have modeled the spectra associated with the atomic, 
molecular and ionized phases of the gas, using different as- 
sumptions for the radiation field, to test the robustness of our 
results. The results obtained by a combination of radiations 
fields along the line of sight confirm the results derived by 
assuming a single radiation field. 

We report a systematic warmer dust temperature in the ion- 
ized phase compared to the atomic and molecular phases. The 
inferred emissivity in the diffuse ionized gas is 2.3 X 10~ 26 
cm 2 /H at 160 fim, which is lower than the Galactic values 
by a factor higher than the metallicity ratio between these 2 
galaxies. This result might suggest different properties of dust 
in the ionized gas of the LMC, compared to our Galaxy. We 
also find, for all cases (1, 2, 3), a significant decrease of the 
PAH relative abundance in the ionized phase compared to the 
other phases, with a larger difference between the phases for 
cases 2 and 3. We interpret this result as due to PAH destruc- 
tion, probably caused by the increased radiation field in the 
ionized phase, although the origin of this phenomenon is still 
under investigation. 

At the same time, the molecular phase appears to favor 
the survival of the PAHs. In addition, when one compares 



bright HII regions with typical HII regions, the ionized phase 
shows an enhancement of the VSG relative abundance by 
more than a factor 2. We also find an important increase of 
the NIR continuum in the ionized phase with respect to the 
other gas phases, which does not seem to correlate directly 
with PAH emission. By comparing results derived using the 
Ha SHASSA map and the Parkes radio data, the estimated 
extinction in Ha for bright HII regions is found equal to 0.5 1 
mag (or A(V)=0.63 mag). 

Finally we observe, for all gas components, a systematic 
increase of the dust emission going from the diffuse medium 
to bright HII regions. 

The origin of the destruction of the PAHs in the ionized 
phase (especially, in HII regions) is still not fully understood. 
Moreover, the depletion of this species is not a systematic 
effect, suggesting that changes in the PAH properties could 
also occur from one HII region to another. A forthcoming 
publication will investigate PAH depletion by means of the 
spectroscopic Spitzer/IRS dat a available in the c ontext of the 
SAGE-Spec Legacy survey (Kemper et al. 2010). This work 
will allow us to study variations of the PAH emission at small 
angular scale within individual HII regions. 
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